Soft repulsive mixtures under gravity: brazil-nut effect, depletion bubbles, boundary 

layering, nonequilibrium shaking 
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Abstract 

A binary mixture of particles interacting via long-ranged repulsive forces is studied in gravity by computer simulation and 
theory. The more repulsive A-particles create a depletion zone of less repulsive 5-particles around them reminiscent to a 
bubble. Applying Archimedes' principle effectively to this bubble, an A-particle can be lifted in a fluid background of B- 
particles. This "depletion bubble" mechanism explains and predicts a brazil nut effect where the heavier A-particles float on 
top of the lighter i^-particles. It also implies an effective attraction of an A-particle towards a hard container bottom wall which 
leads to boundary layering of A-particles. Additionally, we have studied a periodic inversion of gravity causing perpetuous 
mutual penetration of the mixture in a slit geometry. In this nonequilibrium case of time-dependent gravity, the boundary 
layering persists. Our results are based on computer simulations and density functional theory of a two-dimensional binary 
mixture of colloidal repulsive dipoles. The predicted effects also occur for other long-ranged repulsive interactions and in three 
spatial dimensions. They are therefore verifiable in settling experiments on dipolar or charged colloidal mixtures as well as in 
charged granulates and dusty plasmas. 

PACS numbers: 82.70.Dd, 61.20.Ja 



I. INTRODUCTION 

Gravity or centrifugation is commonly used to sort and 
separate different particles out of a mixture^'^ but the un- 
derlying microscopic (i.e. particle resolved) processes of 
mixing and demixing under settling are still debated^i^. 
It is known for long time that shaken or vibrated granular 
mixtures can exhibit the "brazil-nut effect", namely that 
the heavier particles are on top of the lighter ones^i^. The 
details and parameter combinations for the brazil-nut ef- 
fect to occur are still discussed^rJ^. The brazil-nut effect 
may even contradict Archimedes' law, which governs the 
equilibrium density profiles of molecular mixtures and 
colloidal solutions by the buoyancy principle. 

Colloidal mixtures are valuable model systems to ex- 
plore gravity effects on the particle scalei^"— both in 
equilibrium and nonequilibrium. Sediments of binary 
charged mixtures are commonly used to determine the 
phase behavior——. In highly deionized charged col- 
loidal mixtures, the sedimen t^^i^^ can split into sepa- 
rated layers due to counterion lifting, a phenomenon re- 
ferred to as " colloidal brazil nut effect" . The separation 
of binary hard-sphere mixtures was explored using the 
equation of state and separation of the two species was 
predicted in line with Archimedes' principle^^ — . Re- 
cent experimental real-space studies on soft repulsive col- 
loidal mixtures in three dimension^ clearly showed that 
the buoyancy principle is not violated. Colloid-polymer 
mixtures are known to phase-separate under gravity in 
equilibriu m^^i^^ . 

Nonequilibrium studies, on the other hand, include 
the dynamics of the settling process on the particle 



scalei^ii^ after quickly turning the sample upside down, 
the enforcement of crystal growth on a patterned tem- 
plate under gravit y^^i^^ , spatially varying temperature 
fields^^ and novel zone formation in sedimenting colloidal 
mixtures^^. Interestingly, for colloidal mixtures, there 
are only few studies where gravity is changed periodi- 
cally in time^ which may be considered to be the col- 
loidal analogue of granulate shaking. 

In this paper, we consider a binary mixture of particles 
interacting via long-ranged repulsive forces in gravity. 
We use Monte Carlo and Brownian dynamics computer 
simulation and mean-field density functional theory to 
predict the equilibrium density profiles and the nonequi- 
librium response of the system to oscillatory gravity. The 
more repulsive particles are referred to as A-particles 
while the less repulsive particles are the 5-particles. A- 
particles create a depletion zone of small particles around 
them reminiscent to a bubble. Applying Archimedes' 
principle effectively to this bubble, an A-particle can be 
lifted in a fluid background of 5-particles. This "de- 
pletion bubble" mechanism results in a brazil nut effect 
where the heavier A-particles float on top of the lighter 
5-particles. It also implies an effective attraction of an 
A-particle towards a hard container bottom wall which 
produces a boundary layering of the A-particles. If the 
direction of gravity is periodically inverted causing per- 
petuous mutual penetration of the mixture in a slit ge- 
ometry, a similar stable layering emerges as a nonequilib- 
rium phenomenon. We emphasize that our effects do not 
occur for short-ranged interactions like for hard sphere 
mixturesS^. This is only the case when the (effective) 
hard sphere interaction diameter differs largely from that 
of the actual particle size. Therefore, it is the softness of 
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the repulsion which is relevant here. 

Our results are obtained for a two-dimensional bi- 
nary mixture of colloidal repulsive dipoles. Therefore, 
our simulation results can be verified in real-space mi- 
croscopy experiments of two-dimensional superparam- 
agnetic particle^— see also^i^ for alternative set- 
ups. An external magnetic field induces repulsive dipole 
forces. The gravity can either be realized by tilting the 
droplets or by applying a laser light pressure on the 
sample^. The predicted effects also occur for other long- 
ranged repulsive interactions and in three spatial dimen- 
sions. They are therefore verifiable in settling experi- 
ments on dipolar or charged colloidal mixture^ as well 
as in charged granulates^ and dusty plasma a^^i^^ . 

The paper is organized as follows: in section II, we de- 
scribe the model and the simulation technique applied. 
In section III, we discuss the density functional theory. 
Results are presented in sections IV and V and we con- 
clude in section VI. 



II. MODEL AND SIMULATION TECHNIQUE 

The system consists of a suspension of two species of 
point-like super-paramagnetic colloidal particles denoted 
as A and 5, which are confined to a two-dimensional 
planar interface. These particles are characterized by 
different magnetic dipole moments Ma and Mb, where 



M = Mb /Ma 



(1) 



is the dipole-strength ratio. The dipoles are induced by 
an external magnetic field H according to Mi = XiH 
{i = A,B)^ where Xi denotes the magnetic susceptibility. 
The magnetic field is applied perpendicular to the two- 
dimensional interface containing the particles. In the fol- 
lowing, the dipole-strength ratio M is fixed to 0.1, cor- 
responding to recent experimental samples^i^i^. The 
relative composition X = Nb/{Na + Nb) of B particles 
is fixed at 50%; hence we are considering an equimolar 
mixture. The particles are exposed to an external poten- 
tial Vext,i{f^) which is a combination of gravity and the 
hard bottom wall and is given by 



rriigy 
oo 



for ^ > 
otherwise 



(2) 



Here, is the buoyant mass of particle species i {i = 
A^B). Gravity acts along the —y direction. We charac- 
terize the mass ratio by the dimensionsless parameter 



m = mB/TTiA 



(3) 



The particles interact via a repulsive pair potential of two 
parallel dipoles of the form 

Ui^{r) = ^MiM,/r^ = l^^.^jH^r^ = A, B), 

(4) 



where r denotes the distance between two particles in 
the plane. For this inverse power potential, at fixed 
composition X and susceptibility ratio Xb/xa, all static 
quantities depend solely^^ on a dimensionless interaction 
strength (or coupling constant) 



r = 



Mo Xa^^ 
47r kBTR 



(5) 



where /c^T is the thermal energy and Ia = kBT / {rriAg) 
the gravitational length of A particles, which we employ 
as a unit of length. 

For time-dependent gravity ("shaking") we consider 
the external potential 



mig{t)y for < y < Ly 
oo otherwise 



(6) 



which embodies a time-dependent gravity strength g{t) 
in a finite slit of width Ly. This is conveniently modelled 
to be a stepwise constant function of time: 



g{t) 



g for 
-g for 

n = 1,2, 



n 
n 



1 < t/To <n-^ 
^ < t/To < n 



(7) 



introducing a time period Tq and a "swap fraction" < 
'i^ < 1. Note that for the time-dependent gravity we 
confine the system to a slit of vertical width Ly in order 
to keep the external potential bounded from below for all 
times. 

The particle dynamics is assumed to be Brownian. 
Hydrodynamic interactions are neglected. The Brown- 
ian time scale is set by the short-time diffusion constant 
Da of the A-particles. Knowing that this diffusion con- 
stant scales with the inverse of the radius of a particle, 
Db was chosen such that Db/Da = 1.61 correspond- 
ing to the physical diameter ratio of the experimental 
samples^i^i^. 

We perform standard nonequilibrium Brownian dy- 
namics (BD) and Monte Carlo (MC) computer 
simulations^ — in the canonical ensemble where the par- 
ticle numbers Ni {i = A, 5), the temperature T and the 
area A is fixed. In the static (equilibrium) case, where 
the external potential is not time-dependent, the one- 
particle density field has been calculated via MC sim- 
ulations of Na = 600 A-particles and Nb = 600 B- 
particles, which were placed in a finite rectangular box 
Lx X Ly of area A = (30 x 120) (Z^)^- In the nonequilib- 
rium case, where gravity is time-dependent, we obtained 
the one-particle density field by performing BD simula- 
tions of Na = Nb = 300 particles in a rectangular box 
of size A = (15 x 60) (Z^)^- In both cases, the simulation 
box features periodic boundary conditions in x-direction, 
an aspect ratio L^/Ly = 1/4 and a gravitational load 
Na/Lx = 20// A- The coupling constant F is fixed to 
10. A finite time step St = 10~^r was used in the BD 
simulations, where r = 1\/Da- We denote the density 
profiles as pl^ir) {i = A, B) in the static case and Pi(r, t) 
(i = A, B) in the dynamical case. 
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III. DENSITY FUNCTIONAL THEORY 

Within density functional theory (DFT) the grand 
canonical free energy Q{T, /ha, li^B,[PA{r), pB{r)]) de- 
pending on the temperature T and the chemical poten- 
tials jiA-, l-^B is minimized with respect to the local partial 
one-particle densities Pa{t) and psir). This functional 
can be split according tc^^^^, so that in two spatial di- 
mensions we obtain 

n (T, PA, PB, [pA(r), psir)]) = Jid {[pA{r), pB{r)]) + 



^exc([pA(r),pB(r)]) + 



i=A,B 



d^rpi{r) [Fext,i(^) - Pi] , 



(8) 



where the first term is the free energy of an ideal gas 
-Fid = kBT J2 f ^^rpi{r) [\n{Klp,{r)) - l] 



(9) 



i=A,B 



including the (irrelevant) thermal wavelength of parti- 
cles of species i {i = A^B). As already introduced above, 
Kxt,i(^) is the static external potential acting on parti- 
cle species i. The only unknown part is the excess free 
energy functional Texc, resulting from the inter-particle 
interactions. In order to approximate J-'exc, we use a sim- 
ple Onsager functiona l^^i^^ 



dV /dV%(|r-r'|)p,(r)p,(rO, (10) 



(11) 



consisting of the Mayer /-function 

/i,-(r) = l-e-''"^^W, 



with the interaction potential Uij{r) from equation (j4]) 
and the inverse temperature {3 = {kBT)~^. The On- 
sager approximation is valid at low densities reproduc- 
ing the second virial coefficient of the bulk fluid equa- 
tion of state correctly but is expected to break down 
at higher densities. Unfortunately, unlike for hard-core 
interactional^, no alternative approximation working 
at higher densities is known for soft inverse-power-law 
potentials^. However, we expect the general trends to 
be captured by the theory but not the details of molecu- 
lar layering. 

The equilibrium density profiles p^j^ir) are obtained 
from the minimization condition 



5^[pA{r),pB{r)] 



0. 



p*(^)=pr(^) 



(12) 



It is important to note that DFT is typically formulated 
in the grand-canonical ensemble where the chemical po- 
tentials Pa, Pb are fixed instead of the particle numbers, 
while the simulations are performed in the canonical en- 
semble. We have therefore considered the chemical po- 
tentials Pa, Pb as Lagrange multipliers which fix the total 
line density perpendicular to gravity and matched them 



such that this Hue density coincides with that prescribed 
in the simulations. 

If gravity gets time-dependent, there is a dynamical 
generalization of DFT appropriate for Brownian systems 
which can be derived in various ways^^^^ from the ex- 
act Smoluchowski equation via an adiabatic approxima- 
tion. Within this dynamical density functional theory 
(DDFT), the time-dependent density fields obey the gen- 
eralized diffusion equation 



dpijr.t) 
dt 



/3AV- (p^{r,t)\/ 



Sn[pA{r,t),pB{r,t)] 
Spi{r,t) 



(13) 

with a diffusion constant Di corresponding to particle 
species i = A^B. It is important to note that this equa- 
tion conserves the total density, i.e. provided the starting 
density profiles are matched to that in the canonical en- 
semble, the time evolution given by the DDFT equation 
is canonical. Therefore, the results obtained from DDFT 
can directly be compared to our BD simulations. In our 
case of an external potential which depends only on the 
^-coordinate, we consider only density profiles which are 
independent of x. This is justified far away from surface 
freezing^. 



IV. RESULTS IN EQUILIBRIUM 
A. Colloidal brazil-nut effect 

We performed MC simulations for various mass ra- 
tios < m < 1 and dipolar ratios < M < 1. 
Thereby, we choose A as the heavier and stronger cou- 
pled species^. An example for the partial density pro- 
files p^^{r) (i = A^ B) is given in Fig. [1] where at fixed 
M = 0.1 two different mass ratios m = 0.1,0.5 are 
considered. Interestingly, the MC simulation data show 
quite distinct qualitative behavior for these two cases. 
In Fig. [TJi (m = 0.1) the lighter 5-particles are on top 
of the heavier A-particles as expected, while in Fig. [TJd 
(m = 0.5) the behavior is reversed: here, the heavier A- 
particles are on top of the lighter 5-particles. At first 
glance, this opposite trend is counterintuitive. We call it 
- in some analogy to granulate matter - (colloidal) brazil- 
nut effect. 

DFT data are also included in Fig. [TJ In fact, static den- 
sity functional theory can basically reproduce the partial 
density profiles p^j^ir) {i = A, B) though the comparison 
is not quantitative since the functional is approximated 
by a low-density expression. In fact, as compared to the 
simulation data, the DFT results for the layering spac- 
ings are too large but the contact density of A-particles 
at the bottom wall are well reproduced in DFT (see the 
insets shown in Fig. [1]). Of course, one should bear in 
mind that there is no fit parameter involved in the com- 
parison. 

In order to quantify the colloidal brazil-nut effect, we fol- 
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- p;(DFT) 
-- p;^DFT) 

Pb (MC) 
□ (h^,pX)) 
O (hB,P>J) 



no brazil-nut effect 




FIG. 1: (Color online) One-particle density profiles />^^(r), 
i — A^B obtained from DFT and MC simulation. Mean 
sedimentation heights /i^, i — A, B from Eqn. ()14[) are in- 
dicated by highlighted points for MC simulation data, a) 
M 0.1, m 0.1 (no brazil-nut effect), b) M = 0.1, m 0.5 
(brazil-nut effect). The insets enlarge the behavior at small 
wall distances where a strong peak of the A-particles occurs. 



low the criterion proposed in Ref.^^. We define averaged 
heights hi, i = A, B hy taking the first moment of the 
partial density fields as 



hi 



A,B 



(14) 



In Fig. [TJ the location of the corresponding heights are 
indicated by a large symbol. The brazil-nut effect is then 
defined by the condition 



hA < hs , 



(15) 



which means that on average, the heavier A-particles are 
on top of the lighter 5-particles. 

Within the full parameter range 0<m<l,0<M< 
1, the region separating the brazil-nut effect from the or- 
dinary behavior (no brazil-nut effect) is shown in Fig. [2l 
MC computer simulation data for the separation line are 
given by square symbols. These results were obtained 
by systematically scanning the parameter space. The 
brazil-nut effect occurs preferentially for strong dipolar 
asymmetry and is favoured if the two masses do not dif- 
fer much. DFT results for the phase boundary are also 
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FIG. 2: (Color online) Separation line between the occur- 
rence of the colloidal brazil-nut effect and the absence of this 
effect in the parameter space of dipolar asymmetry M and 
mass asymmetry m. Monte Carlo simulation data (contoured 
white squares) , density functional data (full black circles) and 
the transition line implied by the bubble condition fT8)) (con- 
toured circles) are shown. 



included in Fig. [2] and are in good agreement with the 
simulation data predicting the same trends and the same 
slope of the separation line in the M-m parameter space. 



B. Depletion bubble picture 

We now put forward an intuitive picture for the mecha- 
nism behind the colloidal brazil-nut effect which also pro- 
vides a very simple theory for the separation line. This 
picture is based on the observation that when surrounded 
by a fluid of B-particles, a single A-particle creates a cir- 
cular void-like space around it which is free of 5-particles 
due to the strong repulsion between A- and 5-particles, 
see the highlighted area in Fig. [3l This depletion bubble 
is firmly attached to the A-particle. Applying the buoy- 
ancy criterion to the whole bubble, the effective weight 
per area of the A-particles is strongly reduced. Assum- 
ing a homogeneous fluid density pB surrounding a bubble 
of radius R around an A-particle, the brazil-nut condi- 
tion for lifting the A-particle can be stated as a buoyancy 
criterion 

< rriBpB (16) 

We estimate the radius by a low-density argument, 
where the density profile of 5-particles around an A- 
particle fixed at the origin is Pb {—UAB{r)/kBT) = 
exp (— MF/^/r^) such that a reasonable assumption 
for the radius R of the depletion zone is given by the 
separation where the A^-interaction energy equals /c^T, 
i.e. 



R = lA{Mr) 



1/3 



(17) 
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FIG. 3: (Color online) Simulation snapshot depicting A- 
particles (large red spheres) embedded in the fluid of B- 
particles (smaller green spheres) in the bottom region of the 
sample. Notably, A-particles are surrounded by a circular de- 
pletion zone reminiscent of a bubble (area indicated in faint 
red) . Particles are displayed as spheres of finite radii for clar- 
ity only. 




FIG. 4: (Color online) Simulation snapshot for M — 0.1, 
m = 0.5 showing the marked-off bottom layer of heavy A- 
particles (large red spheres) at ^ = beneath the fluid of 
light B-particles (small green spheres). The arrow indicates 
direction of gravity, —y. 



Upon insertion in Eqn. ([T6|) . this yields the brazil-nut 
condition 

m> (n{Mr)hlpB) ^ (18) 

Thereby, a simple estimate for the separation line is pro- 
vided. The only unknown parameter entering in Eqn. 
(p!8|) is the averaged density pB- We have used simula- 
tion data to determine pB as the effective density at a 
distance 2R: 

PB = P'^{2R) (19) 

The resulting separation line is included in Fig. [2j De- 
spite its simplicity, the depletion bubble picture describes 
the simulation data pretty well. Clearly, the depletion 
bubble is induced by the soft long-ranged repulsion and 
is therefore missing for pure hard sphere mixtures where 
neighbouring particles are at contact. However, when 
the interaction is mapped onto a substitute interaction 
core with an effective diameter^^, all the traditional sedi- 
mentation is qualitatively contained in this effective hard 
sphere mixture. 



C. Boundary layering and effective interaction be- 
tween an A-particle and the bottom wall 

We finally discuss the implication of the depletion bub- 
ble on the layering of A-particles close to the hard bottom 
wall of the confining container (at ^ = 0), see again the 
insets of Fig. [TJ The strong layering is clearly demon- 
strated by an actual simulation snapshot shown in Fig. 
m If a single A-particle is fixed at a given distance from 
the bottom wall, its depletion bubble is reduced since 
the void space is cut by the hard wall, see the sketch in 
Fig. [51 Note that the A-particle is point-like so that in 
principle, it can approach the wall very closely. If the 



A-particle is close to the wall, the void space is half of 
the full circle in the bulk (situations I and III in Fig. [5]). 
If the height y of the A-particles increases, the depletion 
bubble area A{y) grows. Assuming a constant depletion 
bubble radius A{y) is given analytically as 



A{y) = R^h~ ^^ccos(|) + 1^1- (i)' foTy<R 
[tt otherwise. 

(20) 

The growing bubble size causes two opposing effects: 
first, in order to increase the depletion bubble area, work 
against the osmotic pressure ps of the fluid 5-particles 
is necessary. Assuming that is constant in the small 
height regime, this work equals pB{A{y) — A{0)) and gives 
rise to an effective attraction of an A-particle close to 
the wall (situation II in Fig. [SJ. In fact, this attrac- 
tion is similar to the depletion attraction in the ordinary 
Asakura-Oosawa-Vrij model of colloid-polymer mixtures 
near a hard wall^i^ although there is a finite (physical) 
colloidal diameter in this model. 

The second effect resulting from the increasing bubble 
size A{y) is a change in the effective buoyancy. The ef- 
fective buoyant force is given by —mAg-\-pB^B9A{y). If 
the bubble containing the A-particle is lighter than the 
surrounding 5-fluid, this term is repulsive with respect 
to the wall and therefore opposed to the depletion force. 

Combining these two effects, we gain the following an- 
alytical expression for the depletion potential V{y) be- 
tween a single A-particle and the wall: 



(pB{A{y)-7rRy2)^mAgy 
V{y) = I -PBTUBg dy'A{y') for ^ > (21) 

[oo for y < 
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FIG. 5: (Color online) Schematic illustration of the 
met ast able trapping of A-particles at ^ = 0, leading to the 
formation of a boundary layer. A-particles are represented by 
large red spheres, while B-particles are depicted as smaller 
green spheres. The solid line indicates y = 0, while the orien- 
tation is analogous to Fig. 2] The depletion zone surrounding 
A-particles is indicated by a dashed outline. 



where the integral can be calculated as 
dy'A{y') =7rR^y - R^y arccos (|) 



,2x1 



(22) 



This expression requires pB as an input parameter. In 
order to evaluate V{y), we have determined pB via a bulk 
reference simulation of a pure 5-system (in the absence 
of gravity) at a prescribed number density by using 
the virial expression^. 

We further checked a posteriori whether the radius R 
of the depletion bubble is consistent with that obtained 
from a radially averaged density profile of a bulk B- 
fluid around a single A-particle fixed at the origin. This 
" renormalized" radius R' can be estimated to be the po- 
sition of the first inflection point in this density profile 
of 5-particles (not shown here). Actually, we find values 
for R^ which are a bit smaller than those given by the 
estimate (pT|) which only works at low pB- 

More rigorously, we can define the effective interaction 
between the wall and an A-particle under the presence 
of the inhomogeneous distribution of 5-particles by a 
potential of mean-force^^'^^. For a given altitude y, the 
effective interaction potential Vef[{y) is given by 



Kff(y) 



Jo 



{FAiy))dy 



(23) 



where {FA{y)) is the canonically averaged total force on 
the A-particle in the presence of the 5-particles (which 
by symmetry points in the ^/-direction). FA{y) also 
contains the trivial direct part —rriAg from gravity. 



Computer simulation data for the effective potential 
Vef^{y) are presented in Fig. [6l Two possibilities occur: 
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FIG. 6: (Color online) Particle- wall depletion potential Kff 
from Eqn. ([23} explored by a single A-particle within a fluid 
of B-particles near the system boundary at = versus 
reduced height y/lA. Panel a) shows the potential barrier 
peak at close distances and includes the theoretically pre- 
dicted depletion potential V{y) from Eqn. (|2T]) for parameters 
M — 0.1, m = 0.5 and two different depletion bubble radii: 
R (black circles), R' (black diamonds), b) and c) show Ven 
for M = 0.1 and M = 0.01, respectively. For clarity, curves 
representing different values of m are shifted relatively. 



if the A-particle is much heavier than the 5-particles, the 
potential is fully attractive by a combination of depletion 
attraction close to the wall and gravity. In the brazil-nut 
case, on the other hand, there are three regimes in Vef[{y): 
an attractive regime close to the wall (corresponding to 
situation I, II in Fig. [5j), followed by a repulsive regime 
(situation III in Fig. [5]) and a subsequent attractive 
regime (situation IV in Fig. [5]). The repulsive regime 
is caused by the lift force also responsible for the brazil- 
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nut effect. 

Fig. [6K includes the theoretical prediction of V{y) from 
Eqn. ([2T]) for one set of parameters. There is very good 
agreement if the renormalized value R' for the depletion 
bubble radius is taken (diamonds) while the agreement 
deteriorates for the low-density-expression R (circles). 
This demonstrates that the analytical expression ([2T]) in- 
corporates the basic physics principles. 

Between the short-ranged wall attraction and the lift 
regime, there is an energetic barrier AE typically of the 
order of several thermal energies ksT. Comparing our 
MC simulation results to approximation (|2T]) with renor- 
malized bubble radius, the energetic barrier height proves 
to be in fair agreement. For M = 0.1, m = 0.5, Eqn. (f2T]) 
predicts AE = 8MbT, while we obtain AE = TMbT 
from our MC simulation. Considering the parameter 
combination M = 0.3, m = 0.9, an energetic barrier 
height AE = Ib.bksT follows from Eqn. (j2ll), while 
MC simulation yields AE = d.SksT. Here, this large 
discrepancy can be attributed to the relatively thin layer 
of 5-particles close to the bottom wall. 

If the energetic barrier is interpreted as a static one, 
an A-particle which is initially trapped in the metastable 
minimum close to the wall needs a huge escape time 
to leave this metastable minimum which scales in an 
Arrhenius-like fashion oc exp (A^/Zc^T). However, as 
pointed out by Vliegenthart et al. in Ref^^, the effec- 
tive interaction can strongly fluctuate such that the en- 
ergetic barrier is not static. In this case, the particle 
will escape much more quickly by waiting for a fluctua- 
tion which reduces the energetic barrier instantaneously 
and accelerates the escape process. Therefore, we have 
computed also the fluctuations of the depletion force 
(AF^) = (F^) - see Fig. [71 Near the potential 

barrier, (AF^) increases significantly. The inset of Fig. 
[7] shows the relative fluctuations. They are indeed of the 
order one near the effective potential maximum. This in- 
dicates that A-particles are exposed to strongly fluctuat- 
ing forces when attempting to escape the boundary layer. 
Therefore, although the depletion potential comprises a 
potential barrier in the order of several ksT^ particle 
transitions from the boundary layer to higher altitudes 
occur at much higher frequencies than expected from the 
static activated Arrhenius expression oc ex.p{AE/kBT). 
In fact, this facilitates the boundary layer sampling in the 
Monte- Carlo simulations of many A-particles to a large 
extent and ensures sufficient equilibration. 
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FIG. 7: (Color online) Mean square fluctuation (AF^) versus 
reduced height v/Ia- Inset: Relative fluctuation (AF^)/(F^) 
versus reduced height u/Ia- The parameters are M = 
0.1, m = 0.5. 



were used in the simulations to obtain statistical averages 
which were all sampled from an interacting bulk system. 
This corresponds to an initial homogeneous density field 
in DDFT. The swap fraction is chosen to be 1/4, i.e. 
we consider the case that the time-average of the gravity 
is non-zero. In particular, we discuss the emergence of a 
steady state upon time-periodic gravity. 

The relaxation of an initially homogeneous (but inter- 
acting) fluid of A- and 5-particles towards its periodic 
steady state can be monitored by observing the instan- 
taneous ensemble- averaged total potential energy £^pot of 
the systerr4^: 

^pot(t) =\ I d\f dV ^ u,j{r-r')p,{r,t)pjy,t) 
+ ^ f d^rVe^tA^,t)p,{r,t) (24) 

i=A,B 



V. RESULTS UNDER TIME-DEPENDENT 
GRAVITY (COLLOIDAL SHAKING) 

Finally, we turn to the nonequilibrium situation of 
time-dependent gravity, see Eqn. [6l which is a sim- 
ple model of colloidal shaking. As far as our methods 
are concerned, we now use Brownian dynamics simula- 
tions appropriate for colloids and dynamical density func- 
tional theory. Various different starting configurations 



This quantity is shown in Fig. [H indicating that only 
few oscillations are needed to get into the steady behav- 
ior. Due to the homogeneous starting configuration the 
energy oscillation amplitude increases with time. DDFT 
describes all trends correctly and also provides good data 
for the potential energies and the associated relaxation 
time. 

The averaged height, as defined by the first moment of 
the density profile (see Eqn. (fT4|) ). can be generalized to 
a dynamical (time-dependent) quantity hi{t)^ {i = A, B) 
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FIG. 8: Total potential energy Epot/{NA + Nb) per particle 
versus reduced time t/r, using time period To = 8t and swap 
fraction 19 = 1/4. The inset shows the short time behavior. 
The parameters are M = 0.1, m = 0.24. 
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FIG. 10: (Color online) A-particle density profile pA versus re- 
duced height v/Ia upon shaking, obtained by DDFT (curves) 
and BD simulations (symbols). The inset depicts the outer 
boundary peak for both methods. Shaking periods To = O.Sr 
and To = 48t are shown for a swap fraction i9 = 1/4. The 
parameters are M = 0.1, m = 0.24. 



(25) 



The time-dependent heights are another indicative pa- 
rameter which probes the dynamical response of the 
whole systen^^i^. The quantity /ia(^) is shown for two 
shaking frequencies a; = 27r/To in Fig. [9l As a result, the 
relaxation time is mainly scaling with the Brownian time 
r but is rather insensitive to the periodicity Tq. DDFT 
reproduces the trend with respect to increasing the peri- 
odicity To but underestimates the actual heights hAit)^ 
consistent with what we found in the equilibrium case. 
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FIG. 9: (Color online) Time-evolution of the mean A-particle 
height hA{t) versus reduced time t/r for different shaking 
periods To and a swap fraction il) — The parameters are 
M = 0.1, m = 0.24. 

Upon shaking, the boundary layer of the A-particles 



persists. We take a snapshot after the relaxation time at 

tn = (n-^)To, n = l,2,.. (26) 

This is just the time at which the direction of gravity 
is inverted. Comparing the results to the density pro- 
files predicted by DDFT, the persistence of the boundary 
layer is indicated by both methods. Fig. [10] depicts den- 
sity profiles obtained by BD simulations and DDFT for 
various shaking periods Tq. Dynamical density functional 
data are in qualitative agreement with Brownian dynam- 
ics computer simulations results but show the same de- 
ficiencies as in the equilibrium case. This indicates that 
the deviations can be solely attributed to the quality of 
the density functional but not to the additional adiabatic 
approximation inherent in any DDFT. Again, as in equi- 
librium, the amplitude of the outermost density peak is 
in good agreement, see the inset of Fig. [TOl 

The mean heights hi {i = A, B) at times tn are shown 
in Fig. [11] versus frequency u = 2t:/Tq. The dissipative 
response leads to a decay of the peak as a function of 
uj. For ^ 0, we recover a quasi-static equilibrium case 
while for oo, the shaking is so fast that the sys- 

tem does not react at all upon this stimulus. It therefore 
approaches the limit of an ordinary equilibrium fluid mix- 
ture confined between two slits in the absence of gravity 
where layering is much less pronounced than in the pres- 
ence of gravity. The crudest model of colloidal shaking is 
a completely overdamped response to a periodic external 
stimulus. In this case, the response amplitude scales with 
the shaking frequency as 1/cj. A fit which involves the 
1/uj scaling is included in Fig. [TT] and provides a good 
description in the high-frequency regime. The DDFT re- 
produces these trends and provides good data for hB{tn) 
while underestimating hA{tn)- Again, we attribute this 
to the low-density approximation of the functional, where 
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the stronger interacting ^d-particles are treated in a more 
approximative way than the 5-particles. 
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FIG. 11: (Color online) Mean heights hi {i — A, B) at times in 
as a function of cj = 27r/To from DFT (squares and diamonds) 
and BD simulations (circles and triangles). The expected scal- 
ing 1/uj is indicated by dashed lines. The parameters are 
M = 0.1, m = 0.24. 



VI. CONCLUSIONS 

In conclusion, we have explored a two-dimensional 
binary mixture of particles interacting via long-ranged 
repulsive forces in gravity by using computer simulation 
and density functional theory theory. The more repul- 
sive A-particles create a depletion zone (void space) of 
less repulsive 5-particles around them reminiscent to 
a bubble. Applying Archimedes' principle effectively 
to this bubble, an A-particle can be lifted in a fluid 
background of 5-particles. This mechanism also work 
when the A-particles are heavier than the 5-particles 
leading to a colloidal brazil nut effect where the heavier 
particles float on top of the lighter particles. Still the 



buoyancy principle is fulfilled if it is effectively applied to 
the mass per bubble volume. This general finding is in 
accordance with the recent experimental sedimentation 
results on colloidal mixtures by Serrano et al^^. 

Within the depletion bubble picture, an effective at- 
traction of a A-particles towards a hard container bottom 
wall is obtained which leads to boundary layering of the 
A-particles. We have also studied a periodic inversion 
of gravity causing perpetuous mutual penetration of the 
mixture in a slit geometry. This non-equilibrium case of 
time- dependent gravity is similar to shaking. Upon shak- 
ing the boundary layering persists. Our results are based 
on Brownian dynamics computer simulations and density 
functional theory. The brazil-nut and boundary layer- 
ing are very general effects, they do also occur for other 
long-ranged repulsive interactions and in three spatial di- 
mensions. They are therefore verifiable in future settling 
experiments on dipolar or charged colloidal mixtures as 
weh as in charged granulates and dusty plasmas. 

Future work should address other interparticle interac- 
tions. Novel effects are expected for a strongly attractive 
cross-interactions leading to mutual mixing of A- and B- 
particles. These interactions are for example realized in 
oppositely charged suspension a^^i^^ . It would be inter- 
esting to check whether a colloidal brazil-nut effect can 
still be observed in this case. Another option for future 
study is to superimpose more external fields (e.g. an ex- 
ternal electric field) to the gravitational field in order to 
controll the response of the system even more^. The 
latter set-up is relevant for electronic in k^^i^^ . 
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